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ABSTRACT 


Data from 29 months of satellite radiation budget measurements, 
taken intermittently over the period 1964 through 1971, are composited 
into mean month, season and annual zonal 1y averaged meridional profiles. 
Individual months, which comprise the 29 month set, were selected as 
representing t!w best available total flux data for compositing into 
large scale statistics for climate studies. A discussion of spatial 
resolution of the measurements along with an error analysis, including 
both the uncertainty and standard error of the mean, are presented. 
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1.0 INTRODUCTION 

A climatology of the net flux of energy exchanged between Planet 
Earth and space has been computed from radiance and irradiance measure- 
ments taken by Earth orbiting satellites. The net flux is derived as 
a difference between total spectral incoming solar flux and the sum of 
separate measurements of reflected shortwave and thermal Infrared 
exitance. This relationship is shown mathematically as: 

Net = Solar In - Reflected - Thermal 

Net, reflected, and infrared flux are frequently referred to as radia- 
tion budget data or measurements in this report. 

Radiation budget measurements are presented in the form of mean 
month, season, and annual zonal profiles. Zonal averaged data are also 
referred to as mean meridional profiles. The terminology zonal averaged, 
is to be interpreted as an average taken over 360 degrees of longitude 
for any given latitude zone. 

Mean zonal radiation budget profiles are presented as climate 
statistics for use in climate studies. The authors believe that the 
profiles for the period 1964 through 1971 are the best statistics 
available to date. Future measurements from the Earth Radiation Budget 
(ERB) experiment on Nimbus 6 and Nimbus G satellites will augment this 
data set. 

There have been a number of requests from persons involved in 
climate research fot such statistics. This report will provide them 
with the statistics and, at the same time, provide a souice of informa- 
tion about satellite radiation budget measurements with appropriate 
references for those desiring additional detailed information. 
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2.0 DATA COLLECTION 


Mean zonal radiation budget profiles are made up from a collection 
of 29 individual monthly sets. This collection does not correspond on 
a one-to-one basis with the collection of an earlier publication by 
Vonder Haar and Ellis (1974), which emphasized maps of data or with the 
mean set of Vonder Haar and Suomi (1971), which did not contain data 
from the early seventies. The collection in this report is shown in 
Table 1. It includes additional data from the ITOS 1 and NOAA 1 satel- 
lites (Flandets and Smith, 1975) during early 1970 and 1971, and ESSA 7 
data (Mac Donald, 1970) in late 1968 and early 1969. Data excluded from 
this report, yet useful for other purposes, ere the TIROS 4 and 7 satel- 
lite measurements in 1962 through 1964 along with 8 months of experi- 
mental satellite measurements. The TIROS satellite could not sample 
poleward of the 63.5 latitudes because of their orbital inclination 
(Bandeen, et al., 1965). Limited on board tape recorder storage left 
data gaps between some ground readout stations. The sampling deficien- 
cies precluded obtaining representative monthly data. 

Measurements from experimental satellites for April through Novem- 
ber 1965 showed large differences between them and the 29 months of 
remaining measurements. Globally averaged albedoes dropped from 28.5 
percent in March 1965 to 19.5 percent in August 1965. Albedo over 
North Africa was in ? neighborhood of 10 percent for the months of 
May, June, July and August 1965. These are extremely low values and 
thus it seems quite reasonable that the data do not represent true ab- 
solute values. Thus, they were not included in the 29 month data zonal 
average data set. 
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Chronological list of earth-orbiting satellites from which the present radiation 
measurements were taken. The approximate local time at which each satellite 
crossed the equator during daylight hours is given in parentheses. EX = Experimental , 
N2= Nimbus 2, N3 = Nimbus 3, E7 = ESSA 7, II = ITOS 1 and NOl = NOAA 1. 
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About the Averages 

Occasionally a near polar zonally averaged albedo was estimated. A 
criterion applied in computing albedo was that if more than 1 watt per 
square meter of incoming solar flux fell into a latitude zone, then 
there should be a reflected flux. Whenever this criterion was not sat- 
isfied, an estimated albedo was assigned to the zone. The criterion 
was not satisfied in a few low i'lumination cases in latitude zones 
bordering the polar night. Here the satellite measured a very small 
signal in the visible light spectrum, a signal not significantly above 
noise in the satellite system. 

Values. which were estimated are shown in Table 2. Estimated al- 

bedoes for these low light cases provide a better input to the net flux 

calculation than an assignment of zero to albedo or reflected flux. 

The most severe case for which an assignment was made is September at 

75 nortn latitude. An assignment of 50 percent albedo gave a 73 watts/ 

2 

m reflected flux. If this estimate is off by ± 5 units out of 100 or 
2 

t 8 watts/m , then it is near the uncertainty in the measurements (dis- 
cussed in Section III). Again this is an extreme case but, there is 
little doubt that an estimated value allows a more representative cal- 
culation of net flux to space than would be obtained by calling a 
missing albedo zero. 

Estimated albedos were not carried through in computing annual 
average albedo. Since measured albedoes are available for many months 
of the data, there is doubt as to whether estimated values would add 
to the representativeness cf annual average albedo. However, estimates 
were included in computing mean season albedo at 65 and 75 south lat- 
itudes in the June- July-August season, and at 85 north latitude in the 
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TABLE 2 


ESTIMATED ALBEDO VALUES 
FUR LOW INSOLATION CONDITIONS 


Month 

Latitude 

A1 bedo 
(Percent) 

Reflected Flux 
Density (Watts/m^) 

June 

-65 

50 

2 

July 

-65 

50 

6 

August 

-65 

bO 

28 


-75 

50 

4 

September 

85 

56 

49 


75 

50 

73 


-75 

60 

50 


-85 

64 

13 

Scdson 




June, July, 
August 

-75 

50 

12 


-65 

50 

1 

September, 

October, 

November 

85 

56 

16 
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September-October-November season. It was necessary to use estimates 
in these seasons since all months in the seasons had missing albedoes 
at such latitudes. 

About the Spatial Resolution 

Resolution of both the measurement and the grid map must be con- 
sidered. Two types of sensor measurements comprise this data set: 
scanning radiometers and wide angle or flat plate disc sensors. 

The scanning radiometers are medium resolution radiometers (MRIR) 
on board Nimbus 2 and 3 satellites. The field of view of the radio- 
meters varies from 50 km of great circle arc distance at nadir to 110 
km at an angle of 40° from nadir (Raschke and Bandeen, 1970). 

All remaining satellite measurements comprising this data set are 
from flat plate disc sensors with a field of view of 180° or 2v sterad- 
ians of solid angle. The solid angle subtended by the Earth at the 
satellite is a function of satellite height only. Thus the spatial 
resolution of a flat plate sensor is dependent on height alone. This 
resolution varies from 53° of great circle arc (5,9U0 km) for lower 
orbiting experimental satellites to 70° (7,770 km) for higher orbiting 
ESSA, ITOS, NOAA satellites. 

If only total power received at the sensor is considered, one may 
be mislead as to measurement resolution of a flat plate sensor. By 
considering a smaller area on the earth's surface contributing to 50 
percent power on the sensor one may get a better estimate of sensor re- 
solution. A great circle arc distance on the Earth's surface contr bitt- 
ing to 50 percent of the power incident on a sensor can be calculated if 
one assumes the Earth atmosphere system to be a homogeneous, isotropic 
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reflector and emitter (Appendix A). Half power resolution in terms 
of great circle arc of the Earth's surface is 11. 5° (1,280 km) for the 
lower orbiting satellite to IQ"^ (2, 130 km) for the higher orbiting 
satellites. Thus, the half power area is only 5 to 10 percent of the 
full power area or approximately 25 percent c the great circle arc. 

All of the data are presented at 10° latitude intervals from 85N 
to 85S in this report. Data from higher resolution scanning radiometers 
were averaged over each 10° latitude zone. If half power resolution of 
the flat plate sensors is considered to be an estimate of sensor measure- 
ment resolution, then it is seen that the experimental data are compatible 
with 10° gridding. However, ESSA, ITOS, and NOAA data, which comprise 
just 10 monhts of our 29 month data set, are much smoother and more re- 
presentative of flux measurements over 20° latitude bands. Users of mean 
statistics presented here should b^; aware of the flux measurement resolu- 
tion. The numbers and graphs should be considered as representing fluxes 
from 10 to 20° latitude zones. 

Flat pla.e data have been reduced from satellite height (hs) to some 
reference height above the Earth's surface (ho). Tiie ho values vary from 
30 km for experimental satellites, 0 km for ESSA 7, and 10 km for ITOS 1 
and NOAA 1; the difference over 0 to 30 km has less than a 1 percent 
effect on the reduced flux value. It must be kept in mind that a reduc- 
tion to some ho is not a deconvolution process which considers inhomo- 
geneously distributed radiation sources in the sensor field of view. 
Instead, homogeniety and isotropy are assumed so that simple geometry 
allow a reduction. It must be noted that a sensor does measure signifi- 
cant anisotropic radiance outside the geometrically reduced field of 
view. The reduction is not too bad when working with time averaged data 
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since transient cloud patterns tend to promote a homogeneous target. 
However, there are certain standing inhomogenities present in time 
average fluxes (primarily due to ice-snow fields, continent-ocean 
distribution, and stationary cloud systems) which preclude simple 
geometric data reduction to an arbitrary reference level . ho. 

3.0 ERROR ANALYSIS OF THE MERIDIONAL PROFILES 

Uncertainty of individual samples can be combined in«*o an un- 
certainty of the mean value. Additionally, standard deviation in net 
flux values can be computed for each monthly time period from idiich a 
standard error of the mean can be calculated. A comparison between com- 
puted uncertainties and error in the mean estimates allows one to draw 
some conclusions concerning natural time variability about the mean. 

Uncertainty in the Mean 

Measurement and data reduction uncertainties are not always well 
known. Each uncertainty is considered qualitatively, at least, as being 
composed of random and systematic errors. The uncertainty due to random 
errors can be minimized by sampling frequently in both space and time. 
Systematic errors, if known in sign and magnitude, can be removed from 
the data. However, some are not known and, therefot'*; r’nnot be removed. 
Individual identifiable, but not necessarily quantitative, uncertainties 
are discussed as follows. 

1) The "solar constant" has been taken as 1360 w/ra^ after Drummond 
et al.(1968). A total uncertainty in the solar constant is es- 
timated to be + 1.5 percent after Thekaekara (1975). 

2) Calibration of sensors and traceability of the calibration to 


primary standards. 
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3) Unaccountable degradation of the sensor in space. 

4) Diurnal sampling bias since all of data are taken from sun 
synchronous satellites, i.e. , satellites which sample at 
the same local sun time each day. Thus, the effects of 
diurnal cloud variations are not measured. 

5) Smoothing in space by flat plate sensors so that a grided 
value represents a measurement for some larger area than 

the grid spacing. This was discussed in the previous section 
on resolution. 

6) Corrections applied to MRIR scanners on Nimbus 2 and 3 to 
account for anisotropic Earth-atmosphere reflections to space. 

7) Parameterization applied to Nin^s 3 longwave spectral 
radiances along with limb darkening parameterization applied 
to both Nimbus 2 and Nimbus 3 MRIR to obtain total longwave 
flux to space. 

8} An assumption of zero net planetary radiation balance applied 

in ESSA 7 data reduction necessary to esolve reflected fluxes 

to space. This assumption becomes less restrictive for longer 

time averaging intervals. Absolute error in net radiation may 

2 

be as large as i 10 watts/m when averaged over a month for 
ESSA 7 data. 

9) All monthly sets have time sampling voids so that a monthly 
mean sample is not quite a true mean. Some monthly samples 
have spatial sampling voids caused by inadequate onboard 
tape recorder storage between satellite ground readout sites. 
Others are due to low signal -to-noise ratio in low light sit- 
uations near the solar terminator on the Earth. 
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10) Natural year-to-year variability of the target which might pre- 
clude the mean of a few monthly samples being a representative 
estimate of a climate mean. 

Conservative estimates of the total uncertainty in incoming solar, 
albedo and infrared exitance are as follow: 

Solar insolation: o = - 1.5 percent 

Albedo: o = t S percent (+ 0.05 x Albedo) 

Infrared Exitance: o = t 5 percent 

Uncertainty in the solar constant of * 1.5 percent is from Thekaekara 
(1975). Uncertainty of t 5 percent in albedo and infrared exitance is 
quite conservative when one considers just uncertainty in sensor cal- 
ibration and degradation which is 2 to 3 percent. However, if we con- 
sider all of the uncertainties in our list, then t 5 percent is not too 
rigid. 

Uncertainty in net radiation has been computed considering effects 

of both dependent and independent errors (Appendix B). Tables 3 and 4 

show the computed uncertainties for mean .r.onths, mean seasons, and mean 

annual net radiation. The very large uncertainties in September are 

due to having just une monthly data set to apply as a mean September. 

The large September uncertainty is not so outstanding in mean season 

uncertainties. Total uncertainty for the mean annual case is less than 

2 

or equal to 10 watts/m at all latitudes. This is not too bad when one 

2 

considers that 1/2 of all the uncertainty or 5 watts/m is equivalent 
to the uncertainty in the global average solar constant. In other words, 
10 watts/m uncertainty in net flux is equivalent to a 3 percent un- 
certainty in the value of the sclar constant if exact Earth flux measure- 
ments could be made. 
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Uncertainty in Mean Monthly Net Radiation (watts/meter^) 
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TABLE 4 

Uncertainty in Mean Seasonal and Mean Annual Net Radiation (watts/meter^) 
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standard Error in the Mean 


Standard error in the mean is an estimate drawn from independent 
random samples of how "good" the mean value is. Standard error in the 
mean is defined as: 

SEM = — 

N 

where a is the standard deviation and N is the number of samples in 
the mean. A sufficient number of monthly radiation budget samples 
are not available to compute SEM on a monthly basis. However, by com- 
bining monthly into seasonal o's, a meaningful statistic can be gener- 
ated. The process of combination is identical to that for combining 
uncertainties (Appendix C). The results in Table 5 show that SEM of 
the polar regions is larger in the fall and winter seasons of each 
hemisphere than is the uncertainty in mean statistics of Table 4. 

Just the opposite, and of lesser extent, is seen in the tropics. The 
large SEM in polar regions indicates that large year-to-year variations 
are in the data in polar regions (70 to 90 latitude), which are larger 
than the uncertainty in the data. They most probably are real inter- 
annual variations. But, one should have less confidence in mean values 
in near polar regions because of large standard error in the mean, 
particularly during the fall and winter seasons. 
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Standard Error of Season Mean and Annual Mean Net Radiation (watts/meter^) 
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4.0 KEAN RADIATION BUDGET STATISTICS 


The statistics are presented so that 12 mean months are followed 

by four mean seasons and a mean annual set. An average value is tab- 

2 

ulated ^or each 10 degree latitude zone in watts/meter except for albedo 
which is In percent. Solar input to the Earth atmosphere system was 
specified with correct Earth-Sun geometry using a solar constant of 
1360 watts/meter^. 

Column headings are defined as follows; 

NET: net radiotive flux exchange with space 

IR: infrared exitance, or thermal longwave flux loss to space 

ALB: Albedo in percent 

ABS: Shortwave or solar flux absorbed in an Earth-atmosphere column. 

REF: reflected and scattered shortwave, or solar flux to space. 

Each mean set of zonal statistics is followed by a graphical pre- 
sentation of zonal average albedo, infrared exitance and net radiation. 


15 



MEAN JANUARY RADIATION BUDGET (watts/meter^) 
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MEAN JANUARY RADIATION BUDGET 
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MEAN FEBRUARY RADIATION BUDGET (WATTS/METER^) 
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APPENDIX A 


Radiant Flux. E^. at satellite height from the Earth Is: 




2it 


m 


L(e,(t>)cose sine de dt 


(Al) 


0 ' 0 


where ♦ and define the right circular cone tangent to the 
earth and subtended by the earth at the satellite and L(6.4>) Is 
radiance. 


If the target Is assumed to be homogeneous and Isotropic, L = L(e,<t>), 
then Al can be Integrated to yield: 

Ej = 2n L sin^ (A2) 

The Interior great circle arc subtended at the earth by this cone Is: 
a'^ = 2o(n = 2(it/2- 6j^) (A3) 

From the geometry in figure Al It Is seen that: 

where h Is the height of the satellite above Earth's surface and 
R is the radius of earth. 


Half power area of the earth, as seen by a flat plate sensor on board 
an earth orbiting satellite at height h is, be equating A2: 
h (2it L sln^ * 2i- L sln^ 0^ 

solving for gives the half angle of this new cone 
0^ * s1n‘^(^s1n Om) 


(A5) 




Figure A1 




From simple geometry, the half power great circle arc at the 
Earth's surface is seen to be: 

= 2(t^ = 2^7t - (y+0^)J (A6) 

where y f>^om the "law of sines" is 
Y = sin ^ sin6|^) 

and the great circle arc distance is R. 


The flux at satellite height (hs) can be adjusted to some other 
reference height (ho). The reference height may be the Earth's sur- 
face (ho = 0) or some height representative of the top of the Earth's 
atmosphere (ho=60km). This adjustment is strictly a geometrical ad- 
justment to the flux with no allowance made for the inhomogeneous and 
anisotropic nature of L(<ji, 0 ) over the target as seen at satellite 
height. 


The geometrical multiplier, 3, using A2 is: 


- / R»tis \ ^ 


so that flux at the new reference level is: 



Table A1 gives various values of the factors for each of the satellites 


with flat plate sensors. 



TABLE Al: GEOMETRY FACTORS FOR SATELLITES WITH FLAT PLATE SENSORS 
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APPENDIX B 


The equation for net radiation is: 

N = (1 - A) S - I (Bl) 

where each term is 
A, albedo 

S, solar flux at the top of the Earth's atmosphere 
I, infrared exitance or thermal radiation. 

The total uncertainty in the net can be expressed by applying a 
Taylor series expansion to small departures from the true value: 



The first three terms on the right-hand side treat the uncertainty of 
each component independently. The covariance or correlated uncertainties 
are in brackets. It is known that uncertainty in the solar constant is, 
to some extent, negatively correlated with uncertainty in albedo since 
albedo is not measured directly, but is computed as function of solar 
insolation. Thus, the first term in brackets is negative. The second 
term in brackets can be eliminated since uncertainties in solar constant 
and infrared exitance are independent. The third term in brackets tends 
to be negative since uncertainties due to sensor degradation in a combined 
system of black and white sensors are negatively correlated. However, not 
all of the data in this report were taken from such a combined system. 

The effect of the two negative terms is to reduce total uncertainty in 
the net. If we consider a worse case when all uncertainties are independent 



we can simplify B2 to an expression cotrmonly referred to as, “the law 
of propagation of independent errors." 



This expression has been applied to each monthly set of data. In 
another form (B3) can be applied to obtain uncertainty in monthly means 
if each monthly set, or each satellite measurement, has an uncertainty 
independent of measurements from other satellites so that 



where n is the number of monthly data sets in each mean month. 

Uncertainty in the seasonal mean must consider that uncertainty in 
each of three mean months of a season are not independent since measure- 
ments from the same satellite frequently appear in each mean month. There 
fore, an equation in the form of B2 must consider a.l six terms with posi- 
tively correlated uncertainties. Such an equation can be simplified to be 



where K is the number of mean months in a mean season (k < 3). 

An equation identical to B5 can be applied to each mean season un- 
certainty to obtain uncertainty in the annual mean net radiation. 



APPENDIX C 


Standard Error In the Mean 

A mathematical formulation for computing standard error in the mean 
(SEM) is developed from dependent and independent error analysis as pre- 
sented in APPENDIX B. 

The SEM of the net radiation for each mean month considering in- 
dependent data sets is: 




error in the mean for four seasons and the annual case. 



